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I. INTRODUCTION 



Currently the observational evidence coming from su- 
pernovae studies cosmic background radiation fluc- 
tuations Q, and baryon acoustic oscillations 0], set a 
strong case for a cosmological (concordance) universe 
model composed by nearly 70 percent of a mysterious 
component called dark energy, responsible for the current 
accelerated expansion, nearly 25 percent of dark matter, 
which populate the galaxies halos and a small percent 
(around 4) composed by baryonic matter. The nature 
of these two dark component remains so far obscure [H • 
In the case of dark matter, we have a set of candidates 
to be probed with observations and detections in parti- 
cle accelerators, however the case of dark energy is more 
elusive. This component not only has to fill the universe 
at the largest scale homogeneously, and to have the ap- 
propriate order of magnitude to be comparable with that 
of dark matter, but also has to have a negative pressure 
equation of state, something that was assumed first to 
explain inflation in the early universe, inspired by high 
energy theories of particle physics, but which seems to be 
awkward to appeal for, at these low energy scales. The 
alternative way used to account for the cosmic accelera- 
tion is to consider a modification of general relativity at 
large scales Q. However it seems to be a difficult route 
to follow, considering the subtleties in the interpretation 
of the observations even in the simplest model based on 
general relativity. 

To describe the universe at the largest scale, as is the 
case in cosmology, we have to make certain assumptions 
regarding the local distribution of matter and its char- 
acteristics. As is well known (and explained in several 
classic books @, 0, @|) the energy-momentum tensor 
used in the Einstein's equations is imported from con- 
siderations made in special relativity, through a covari- 
ant generalization supported by the equivalence princi- 
ple. The assumption about homogeneity and isotropy 
leads to consider the matter content of the universe as a 
fluid. In the case of a perfect fluid, defined as a mechan- 
ical medium incapable of exerting transverse stresses, it 
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takes the standard diagonal form 

T"" = {p + p)U^U v -g^p, (1) 

where p is identified as the energy density of the fluid, p 
is the proper hydrostatic pressure, J7 M = dx^/ds are the 
components of the macroscopic velocity of the fluid with 
respect to the actual coordinate system and g^ u is the 
metric of it. The expression |T]) is what is measured by 
an observer at rest in the fluid, who examines an element 
of the fluid small enough so that gravitational curvature 
can be neglected. The conclusion of applying the energy 
conservation relation to this tensor 

= o, (2) 

is that perfect fluids behave adiabatically when examined 
by a local observer, satisfying the equation 

d(pV)+pdV = 0, (3) 

where V is the physical volume of the region considered. 
Is usual to take as the volume V = c?y, where df is a 
physical distance of the order of the region we are inter- 
ested in. Using that df — d c R(t), leads to the well know 
result 

P + 3|(p + p)-0, (4) 

To get this result, it was assumed also that particle num- 
ber is conserved Nft = 0, where = nU^ with n the 
particle number density. 

A way to modify ([3]) without spoiling adiabatic evo- 
lution is by adding a work term due to the change of 
the particle number N. In Q, the authors studied this 
case as an alternative cosmological model, considering 
the change of the particle number, assuming that such a 
correction have to be considered during matter creation. 
This work actually suggested for the first time, a way to 
incorporate the particle creation process in the context of 
cosmology, in a self consistent way. In fact, the original 
claim made by Zeldovich (loj . that gravitational parti- 
cle production can be described phenomcnologically by 
a negative pressure, is realize here in a beautiful wa y. A 
covariant formulation of the model was presented in jllT ] . 

The context where these considerations would be im- 
portant were mentioned to be; the steady-state cosmo- 
logical model, the warm inflationary scenario and within 
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the standard inflationary scenario, during the reheating 
phase. Based on this work, the studies of cosmological 
models with matter creation [T3] were initiated, which 
were rapidly recognized as be potentially important to 
explain dark energy [13] . In particular, in Ref. [1J| the au- 
thors established a model where dark energy can be mim- 
icked by self-interactions of the dark matter substratum. 
Within the same framework, models of interacting dark 
energy and dark matter were proposed . Actually, in 
these studies is possible to have consistently, a universe 
where matter creation proceeds within an adiabatic evo- 
lution. More recently, Lima et al.,[lj| have presented 
a study of a flat cosmological model where a transition 
from decelerated to accelerated phase exist. They ex- 
plicitly show that previous models considered [13], does 
not exhibit the transition, and study the observational 
constraints on the model parameters. 

In this work I consider non flat models where by mod- 
ifying the first law, enable us to explain the current ac- 
celeration of the universe expansion through a fictitious 
pressure component coming from changes in the total 
particle number. In this way, the dark energy compo- 
nent is basically an effect due to the dark matter creation 
process in the universe, fact that enable us to explain eas- 
ily the cosmic coincidence problem; it is not strange to 
have a similar contribution from these two component, 
because they are just one component; dark matter. I 
also demonstrate that a non constant matter creation 
rate can indeed lead to a transition from a decelerated to 
an accelerated regime even in the case of non flat geome- 
tries. In the next section I derive the equations of motion 
in the case of matter creation. Then, I open section III 
with a very simple model with a non constant creation 
rate, that resemble the ACDM model. After this, I dis- 
cussed the case of the transition from a decelerated to an 
accelerated expansion in non flat models. 



II. MATTER CREATION 

The incorporation of the variable N into the cosmo- 
logical scenario, modifies the first law as 



d(pV) +pdV - (h/n)d(nV) = 0, 



(5) 



where h is the enthalpy (per unit volume), n = N/V is 
the number density and p is the energy density. The im- 
portant thing to stress here is that the universe evolution 
continues to be adiabatic, in the sense that the entropy 
per particle remains unchanged ((7 = 0). The extra con- 
tribution can be interpreted as a non thermal pressure 
defined as 



Pc 



h d{nV) 
n dV ' 



(6) 



This is the source that produces the current acceleration 
of the universe expansion. Once the particle number in- 
creases with the volume (dN/dV > 0), we obtain a neg- 
ative pressure. 



Now, because we are considering matter creation, we 
have to impose the second law constraint as 



dS = -d{nV) > 0, 



(7) 



where s = S/V is the entropy density. From ([5]) we find 
that the new set of Einstein equations are 



8ttG 



1 r+ & = 



p= ~{p + p) 
n 



(8) 



(9) 



and from Q we have to supply an extra relation between 
n and the Hubble parameter 



3Hn = 3H(3, 



(10) 



where (3 is a definite positive function. The set of equa- 
tions (|S]), © and (JTUJ) completely specified the system 
evolution. The standard adiabatic evolution is easily re- 
covered: setting (3 = implies that h/n — —3H, which 
leads to the usual conservation equation from ((9|). It is 
interesting to note that Eq. ([9]) enable us to determine the 
pressure; for example if p — win Eq.f9|) implies p = 0, and 
furthermore if p — aT A and n — bT 3 implies p = p/3. 



III. MODELS 



Let us assume first here that 

dN -6 



(11) 



The arguments in favor of this hypothesis would be given 
in short. Let us study here its consequences. If we con- 
sider V — R 3 , relation (JTTJ leads to /3 = /3 . From (JlOj) 
we finds the solution for the density number to be 



n(t) = A) - 



_Cj_ 

R(ty 



(12) 



where C\ is a constant of integration. Assuming p = mn 
as for non-relativistic matter (that leads to que usual 
equation of state p = as we discussed that at the end 
of the last section), we get for the energy density 



p(R) = mPo 



(13) 



that resembles the combined contribution of a cosmo- 
logical constant and dust. The level of fine tuning here 
to obtain an accelerated expansion is relatively smaller 
than in the usual ACDM model, because in this case, 
both terms come from the same function, and we know 
that we can obtain an accelerated expansion phase after 
some time from ©. Note however that the equation of 
state has not been modified, p is still the non-relativistic 
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contribution. We do not have to introduce any exotic 
component - with a negative pressure - to describe the 
current expansion acceleration. Given the current status 
of the dark matter and dark energy problem Q , we can 
use this idea as a possible way to understand it. 

In terms of the exotic dark energy fluid with energy 
density px and pressure px, the solution (TIB")) can be 
expressed as 



with 



= ni JT 



(17) 



where 7 is a constant parameter (is exactly the same let- 
ter introduced in [Hj]). In what follows, a non relativistic 
matter equation of state is assumed (p — 0). 

Introducing p7)l in (fTQ|) leads to the following solution 



px = m(3, 



px = ujpx , 



(14) 



where now there must be two copies of Eq.© with the 
one corresponding to dark energy being 



Px = ~{px +px), 
n 



(15) 



and n stand only for the second term in Eq. (|13p . Of 
course, from Eqs. (JSJ), (fl3|) and (fT5|) we get u = — 1; 
e. L, the ansatz (|TT|) has leads us to the case of a 
purely cosmological constant contribution. In general, 
we can always separate the contribution for "purely" non- 
relativistic matter p m (oc R~ 3 ) from the dark energy. In 
this case p = p m + px, and the pressure can be written 
as 



n(R) = n a 



R 



(18) 



Is evident the meaning of the subscript zero. Because p = 
mn for non-relativistic matter, the Friedman equation 
can be written as 



R 2 R + kR = Ke 3ffo7( *- to) , 



(19) 



where k = 8ttGpoRq/3. Evaluating this expression for 
the reference time t — tg we obtain that 



k = H$R 2 (n - 1), 



(20) 



where we have defined flo = p/p c , with p c — 3Hq/8itG 
being the critical energy density. Using this in (|19[) we 
obtain 



Px 



-m(3 . 



(16) 



We can notice that px is exactly the non thermal pressure 
computed in ([6]) (recall that p m = 0). 

I have to stress here that the main result derived in this 
section means that we have a single contribution, which 
satisfy the non-relativistic matter equation of state, that 
describe both dark matter and dark energy simultane- 
ously, and in this way solving automatically the coin- 
cidence problem. This dark unification mechanism does 
not have the problem studied in [18[ , because in the cases 
studied in that paper, for example the Chaplygin gas 
model [H], the sound velocity of the dark matter is not 
zero, leading to instabilities. That happens because the 
Chaplygin gas interpolates between the equation of states 
for dark matter and dark energy. This does not happens 
here, because there is just one equation of state; that of 
non-relativistic matter. 



IV. DECELERATED/ ACCELERATED 
TRANSITION 

As is well known from observations of supernovae, a 
transition from a decelerated to an accelerated expansion 
occurs in the recent history of the universe. Depending 
on what is used to model dark energy, different redshift 
have been obtained between 0.5 to 1. 

This topics was recently discussed in [l6[ where the 
authors specialized in a flat cosmology, where a explicit 
transition can be achieved from a decelerated to an accel- 
erated expansion. In this section we generalize this work 
to non flat universes, by using a model (see equation (|10|) ) 
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(21) 



Clearly, for 7 = this expression reduces to the standard 
result. 

In order to compare with observations, we shall com- 
pute the deceleration parameter, q = — RR/R 2 , in terms 
of the redshift. By differentiating (|2"Tj) respect to cosmic 
time we find 



(3iJ 7 - H)e 3Ho '<( t - to ') 
'2H[l-n + ft (-7f )e 3ff o7(*-to)] 



(22) 



So, from ([21]) we find the scale factor in term of the cosmic 
time, R(t), and using that R — Rq(1 + z)^ 1 , we can use 
(|2"2"|) to write down q(z). As an example, for CIq = 1.2 
and 7 = 0.4 the deceleration parameter computed from 
(|22p in terms of the redshift is shown in Figure 1. For a 
fixed flo, a crossing redshift exist if 7 > 1/3. Otherwise, 
the universe never made the transition. For Qq < 1, 
the situation is qualitatively similar, in this case with a 
smaller slope. 

In the case of a flat universe, tto = 1 we obtain 



So 



37^) 



(23) 



which coincides with the result obtained in [ll| . 

To complete the analysis, I will briefly discuss the case 
of the interaction term T = 3[3H, discussed also in [l2l ]. 
In this case, the deceleration parameter can be explicitly 
written in terms of the redshift as 



?(*) = 



1-3/3 



n (i + zy-w + 1 - n 



(24) 
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FIG. 1: q(z) versus z. An example of Eg. (|22[) in terms of the 
redshift z. In this case f2o = 1-2 and 7 = 0.4. See the text for 
comments. 
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FIG. 2: Typical profile of the scale factor evolution in the case 
of a oscillatory distribution of matter. Although the global 
behavior is oc R~ 3 , the oscillatory nature of the distribution 
of matter produces wiggles that can be associated to short 
periods of accelerated expansions. The scales are arbitrary. 

If we specialize to a flat universe, Qo = 1, q becomes a 
constant, showing no transition between decelerated to 
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an accelerated regime. However, the non flat case, is not 
different. Actually, the z dependence make the values of 
q to vary with redshift, but there is no redshift for which 
3 = 0. " 

So the challenge would be to estimate the function 
dN/dV enabling us to give the best fit to the obser- 
vations. The first model considered in the last section 
shows us a very simple scenery to follow. Let me pro- 
ceed further and complicate a little bit the model setup. 
For example, we can study the case of a distribution of 
matter which is oscillatory with volume. It means that 
instead of consider (fTTj) . we use dN/dV — f3 cos(V/V c ). 
This implies that matter distribution is characterized by 
a volume V c . In this case, the energy density behaves 
as p = d/R 3 + {(3V c /R 3 )sm{R 3 /V c ), so when the vol- 
ume considered is small enough R 3 -c V c , this solution 
approaches (fllij) . A typical profile of the scale factor evo- 
lution is shown in Fig. 2. 

In this letter I have demonstrated that considering 
small changes in the total number of particles in our 
universe, offer a possible way to understand the current 
accelerated expansion measurements, without using any 
exotic energy component. Also, this scenario explain also 
the cosmic coincidence problem, basically showing that 
these two components, dark matter and dark energy, are 
the of the same nature, but their act at different scales. 
This way of understand the SNIa observations, implies 
that cosmology have a new window to explore the uni- 
verse considering matter creation, encoded in the func- 
tion dN/dV. Although I have discussed some examples 
for this function, neither has been obtained from a fun- 
damental theoretical basis. 
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